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Abstract 

Lie symmetry method is appUed to find analytic solutions of initial-boundary- 
value problems of transient conduction in semi-infinite solid with constant surface 
temperature or constant heat flux condition. The solutions are obtained in a 
manner highlighting the systematic procedure of extending the symmetry method 
for a PDE to investigate BVPs of the PDE. A comparative analysis of numerical 
and closed form solutions is carried out for a physical problem of heat conduction 
in a semi-infinite solid bar made of AISI 304 stainless steel. 

1 Introduction 

Lie symmetry method is a powerful general technique for analyzing non-linear PDEs 
and can be efficiently employed to study problems having implicit or explicit symme- 
tries. Thus it provides most widely applicable technique to find the closed form solu- 
tions of differential equations and contains, as particular case cf. [20], many efficient 
methods for solving differential equations like separation of variables, traveling wave 
solutions, self-similar solutions and exponential self-similar solutions. Since the modern 
treatment of the classical Lie symmetry theory by Ovsiannikov [19], the theory of sym- 
metries of differential equations has been studied intensely and has substantially grown. 
A large amount of literature about the classical Lie symmetry theory, its applications 
and its extensions is available, e.g. [1, 2, 5, 6, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 21]. 
Most of the engineering and physical problems require the PDE to be solved subject 

to suitable initial or boundary conditions. Although there have been some notable 
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contributions in the applications of symmetry method to boundary value problems 
(EVPs) c.f. [1, 3, 7, 8, 19], in general the Lie symmetry method has not been utilized 
in a great deal in obtaining solutions of EVPs of physical significance. One reason 
could be the natural restrictions imposed because of the requirement of invariance of 
initial-boundary conditions in addition to requirement of invariance of PDE under the 
symmetry. This work is concerned with application of Lie symmetry method to obtain 
analytic solution of two standard initial-boundary-value problems of heat conduction, 
in a systematic manner that highlights the systematic procedure of extending the 
symmetry method for a PDE to investigate EVPs of the PDE. 

A description of the initial-boundary value problems, under study, is provided in 
Section 2. Section 3 is divided in two parts, giving an introduction to the solution- 
method followed by its application to obtain exact solutions of the initial-boundary 
value problems described in Section 2. The comparative analysis of numerical and 
analytic solutions for the heat conduction in sohd bar made of AISI 304 stainless steel 
is presented in Section 4. 



2 Description of the initial-boundary value prob- 



We consider test problems related to transfer of heat by conduction. The analysis 
of such problems is required in many physical engineering problems, for example, the 
cooling of electronic equipment, the design of thermal-fluid systems, and the material 
and manufacturing processes. In practice, a major objective of the solution of such 
problems is to determine the temperature field in a medium as a result of either thermal 
condition applied to the boundary of the medium or heat generation within the medium. 
Once the exact temperature distribution is known, the heat flux at any point in the 
system, including the boundaries, can be computed form the Fourier's law. 

The problems studied here are for the transient conduction in semi-infinite solid 
either with constant surface temperature or with constant surface heat flux conditions. 
Precisely, we investigate the following initial-boundary value problems. 



Transient conduction in a semi-infinite solid with constant surface tem- 



lems 



IBVP-l 



perature 



dT 



(2.1) 




with initial and boundary conditions 



\t=o 



Til ^|a;=0 — T\r^^^ — Ti (2-2) 



IBVP-2 

Transient conduction in a semi- infinite solid with constant surface heat 
flux 

with initial and boundary conditions 

dT 

r|t=o = o, -fc— |,=o = go", ru_oo = o (2.4) 

3 Lie symmetry solution of the initial-boundary value 
problems 

The general procedure of applying Lie symmetry method to study IBVPs of PDEs 
requires [4] determination of a one parameter Lie group of transformations that leaves 
the problem invariant, and the utilization of these transformations either to construct 
the invariant solution or to obtain similarity reductions. It can be explained by the 
following steps: 

• Lie symmetries of the PDE 

Determining the symmetry algebra of the governing PDE. 

• Invariance of the boundaries 

Taking the most general symmetry operator X obtained in step 1 and finding 
the conditions under which it leaves the boundaries invariant. 

• Invariance of boundary conditions on the boundary 

Finding the restrictions on X that are imposed due to invariance of boundary 
conditions on the boundary. 

These steps will determine the symmetry operator that leaves the IB VP invariant. 

• Construction of similarity solution or reductions 

Utilizing the similarity variables of the symmetry operator of the IBVP to find 

the similarity reductions and similarity solution of IBVP. 
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Further details about these steps are illustrated in the subsections below. It should 
be noted that the more general the symmetry operator (leaving the IBVP invariant) 
found the more likely it is to lead to the solution of the problem [16]. 



3.1 Lie symmetries of PDE (2.1) 

The method of obtaining the classical Lie symmetries of a PDE is standard which is 
described in detail in many books, e.g. [6, 11, 15, 18, 21]. For instance, to obtain 
Lie symmetries of an equation like (2.1), one considers the 1-parameter Lie group of 
infinitesimal transformations in {x, t, T) given by 

X* = a; + e^(a;,t,T) + 0(e2) 

t* ^t + er{x,t,T) +0{e^) 

T* = T + e(i){x,t,T) + 0{e^) 

where e is the group parameter, hence the corresponding generator of the Lie algebra 
is of the form 

d d 

If X'^l denotes the second prolongation of X then using the invariance condition 

XPl (Ti-ar,,)|T,=„T.. = 

yields an overdetermined system of hnear PDEs in ^, r and called determining equa- 
tions. The general solution of this system determines the generator of the symmetry 
algebra. For most type of problems. Computer Algebra Systems like Macsyma, Maple, 
Mathematica, MuPAD and Reduce can be used either to perform above steps sepa- 
rately or to find the symmetry algebra directly. A detailed survey of software that can 
be used for symmetry analysis is given by Hereman in [14, pages 367-413]. 

The symmetry algebra of PDE (2.1) is well known (in fact it was found by Lie) and 
is spanned by the vector fields 

d d d d did 
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The corresponding optimal system of 1-dimensional subalgebras is given by the follow- 
ing vector fields, [18]. 



^1 + cXg, Xi + X4, Xi — X4, Xi + X5 + cXg, X3 + cXq 

and representatives of the form X^o- This means any invariant solution of PDE (2.1) 
can be found via a suitable transformation of the invariant solutions obtained from 
the symmetry operators in the optimal system. Instead, for our purpose we directly 
find the symmetry operator that preserves the boundary and boundary conditions and 
hence directly leads to the solutions of initial-boundary-value problems (see details 
below) . 

3.2 Solution to IBVP-1 

We consider the general symmetry operator 

X = kiXi + k2X2 + ksXa + k^X^ + k^X^ + k^X^ 

of PDE (2.1) and search for the operator that preserves the boundary and the boundary 
conditions (2.2). 

The invariance of the boundaries x = 0, i = or equivalently 

[X(a;-O)],=o = 0, 
[X{t-Q)]t=o = 



implies 

Hence, X must be 



k^^k2 = ki^Q. (3.1) 



X = k^X^ + + k^X^. 



In addition to the restrictions imposed by Equation (3.1), the invariance of initial and 
boundary conditions i.e. 

[X{T-T,)]t=o = 0, on T = Ti 
[X{T-T,)l=o = 0, onT^T, 

implies we must have 

/cg = = ke. (3.2) 
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Hence the IBVP (2.1), (2.2) is invariant under the symmetry 

^ ^ d d 
X^2t- + x—, 
at ox 

where we have chosen — 1. 

The invariant solution of the problem is constructed by utilizing the transformations 

through similarity variables for X. Solving the characteristic system for XI — gives 

x^ d d 

h — — and I2 — T ets the differential invariants oi X — 2t— + x—. Hence, the 

t at ox 

similarity variables for X are 

2 

e(x, t)^j and V{0 = T. (3.3) 

Substitution of similarity variables in Equation (2.1) implies that the corresponding 
similarity solution of PDE (2.1) is of the form T = V{^) where V{^) satisfies the ODE 

The above equation can be integrated, using the substitution 

di' 

to obtain 



V{0 = ci / -^d^ + C2. (3.5) 



Making the change of variable 
in the above solution yields 



V = 4:Cl^/a J e '"^ dy + C2 
= 2ci ^/a^ • erf (y) + C2 , 



where erf denotes the error function. Hence, from Equations (3.3), (3.5) and (3.6), the 

d d 

exact solution of PDE (2.1) that is invariant under X = 2t— + x— is 

at Ox 



r(x,t) =2ciVav^-erf f^^^^^ +C2. (3.7) 



Imposing the initial and boundary conditions determines 

T —T 
ci = — and C2 = T^, 
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giving the solution 



T{x,t)^{Ti-T,)erf 



X 



2 Vat 



of the IBVP (2.1), (2.2). 

3.3 Solution to IBVP-2 

Taking the general symmetry operator 

X = kiXi + k2X2 + hX^ + kiXi + /CgXg + k^XQ 

of PDE (2.3), we determine the operator that leaves the boundary and the boundary 
conditions (2.4) invariant. 

The invariance of the boundaries x — — Q or equivalently 

[X(x-O)],=o = 0, 
[X{t-Q)]t=o = 

ki — k2 — ki — 0, 

X = k^X^ + hXr, + k^XQ. 

[X(T-0)]i=o = 



yields 



i.e. X must be 



Since 



(3.8) 



implies 



'^'-x^ + fee 1 T = 0, 



a 



(3.9) 



(3.10) 



the invariance of the condition T|f=o = at T = does not impose any restrictions on 
X. 

The invariance of the boundary condition 

dT 
-k— = go 
ox 



x=0 



requires, see [4, Chapter 4] for details. 



ox 



x=0 



dT 

on - k— = go" 
ox 



where X'^l denotes the first prolongation oi X. It follows that 



dT' 

[-6ht + h-k3)— =0 

J x=0 

on 

dT 

+ qo" = 0, 

ox 

hence we must have 

fcs = and ^3 = kg. 

Choosing k^ — ke = 1 provides the symmetry 

d d d 
X ^Xs + Xe^x— + 2t— + T— 
ox at oT 

that leaves the IB VP (2.3), (2.4) invariant. 

To find the similarity transformations that will lead to the solution, the character- 
istics system for XI = is solved. This provides the similarity variables 

e(a;,t) = ^andV(0 = -. (3.11) 

Substituting similarity variables in PDE (2.3) implies that the corresponding similarity 
solution is of the form 

T = xViO (3.12) 

where V{^) satisfies the ODE 

Following a procedure similar to the solution of ODE (3.4) and using Equations (3.11), 
(3.12) we obtain 

T(x, t) = -2ciVie^ - cix^l^ ■ erf (^^1^^ + C2X. (3.14) 

\ a \ 2^0L J 

Imposing the boundary conditions determines 



go" fa , ?o" 



giving the solution 



T[x, t) = 2— W — e 4at + —X < erf 



k \ n k I \2\/ai 

of the IBVP (2.3), (2.4). 



Remcirks 



1. Although the above procedure illustrates that the symmetry algebras admitted by 
boundary value problems, in general, may not be rich enough to obtain symmetry 
solutions or reductions of the BVP, it is worth noting that sometimes the invariant 
solutions obtained from PDE can satisfy boundary conditions that are not left 
invariant by the operator which generated the solution. Examples can easily be 
constructed from above boundary value problems. For instance, we can see from 
above that the invariant solution 

r(x,t)=r,erf 



obtained from the generator 



satisfies the IBVP 



X3 = 2t- + x— , 
at ox 



dT _ d'^T 
dt dx"^ 



with initial and boundary conditions 

T|t=o = ^|a;=2VS 

but the boundary condition 

T = T^erf 
is not preserved by the generator 



1 

^7^ 



1 



x=2^/a 



2t— — 
dt dx' 



2. In case of non-linear problems it is likely that reduced ODE, obtained through 
similarity variables, may not be integrable in terms of known functions. In such 
cases, the reduced BVP of ODE either can be solved numerically or a Lie plane 
analysis, cf. [8], of the reduced ODE has to be carried out to understand the 
solution. 
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4 Comparison of numerical and analytic results for 
a physical problem 

To compare the exact solution of the two boundary value problems discussed in the 
previous sections with numerical results using finite element method, transient heat 
conduction in a semi-infinite solid bar was solved. The solid is made of AlSl 304 
stainless steel and initially at a temperature Tj. The material properties for AISI 304 
steel used in this work are given below: 

Thermal conductivity (k) = 18.2 W/m 
Specific heat (c) = 0.536 kJ/kg 
Density (p) = 7822 kg/m^ 

Thermal diffusivity {a = k/pc) = 4.ME - 3 rri^ / s 



The spatial temperature distribution at different times was calculated when the 
solid is subjected to the following boundary conditions: 



For BVP-1 


For BVP-2 


Initial temperature Tj = 300 '^A' 

Surface temperature Tg = 900 '^K at a; = 


Initial temperature Ti = 
Heat fiux = 5 kW/m^ at a; = 



In order to numerically solve this problem, a transient thermal analysis is performed 
using FEA software ANSYS. The given structure is modeled using 3-D Conduction Bar 
Elements (LINK33). L1NK33 is a uniaxial element with the ability to conduct heat 
between its nodes. The element has a single degree of freedom, temperature, at each 
node point. The conducting bar is applicable to a steady-state or transient thermal 
analysis. A refined uniform mesh is used to model the nonlinear thermal gradient 
through the solid. The length of the model is taken as L = 2m for BVP-1 and 
L = 10m for BVP-2 assuming that no significant temperature change occurs at the 
interior end point during the time period of interest. This assumption is validated by 
the temperature of node at x = L at the end of the transient analysis. 

In the figures given at the end. Figure 1 shows the spatial distribution of tempera- 
ture using closed form solution of the BVP-1. The comparison of this solution with the 
numerical results using FEA at three different times is shown in Figure 2. Similarly, the 
temperature distribution for the BVP-2 from the exact solution is given in Figure-3. 
Figure-4 indicates very good agreement between closed form solution and numerical 
results for the BVP-2. 
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5 conclusion 



This paper presents a systematic approach for applying Lie symmetry method to inves- 
tigate boundary value problems of partial differential equations. Equations for transient 
conduction in semi-infinite solid are considered subject to either constant surface tem- 
perature condition or constant heat fiux condition. Those symmetries are determined 
that leave the partial differential equation, boundaries and the boundary conditions in- 
variant. The corresponding similarity transformations are then utilized to transfer the 
problem to boundary value problem of ordinary differential equations which are solved 
to obtain explicit analytic solutions in both cases. Finally, the results are applied to 
practical problems of heat conduction in a solid bar made up of AISI 304 stainless 
steel. For these particular problems, a comparative analysis of numerical and closed 
form solutions is also carried out. The recovery of well-known solutions of important 
transient conduction BVPs of engineering shows that a systematic application of Lie 
symmetry method to BVPs is an applicable technique and demonstrates its potential 
for investigating BVPs related to practical applications in different fields of applied 
sciences and technology. 
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Analytical SDlution of BVP # 1 




Figure 1: Temperature distribution at different times using analytical solution 
BVP-1 
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Figure 2: Comparison of analytical and numerical results for BVP-1 
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Figure 3: Temperature distribution at different times using analytical solution 
BVP-2 
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Figure 4: Comparison of analytical and numerical results for BVP-2 
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